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^ Abstract 

o 

(N 

Q In this paper, we propose clipping noise cancellation scheme using compressed sensing (CS) for 

orthogonal frequency division multiplexing (OFDM) systems. In the proposed scheme, only the data 
tones with high reliabiUty are exploited in reconstructing the clipping noise instead of the whole data 
tones. For reconstructing the clipping noise using a fraction of the data tones at the receiver, the CS 
technique is applied. The proposed scheme is also applicable to interleaved orthogonal frequency division 
multiple access (OFDMA) systems due to the decomposition of fast Fourier transform (FFT) structure. 

O 

I— ' Numerical analysis shows that the proposed scheme performs well for clipping noise cancellation of 

CnI both OFDM and OFDMA systems. 

> 

^ Index Terms 

Clipping noise, compressed sensing (CS), fast Fourier transform (FFT), orthogonal frequency divi- 
sion multiplexing (OFDM), orthogonal frequency division multiple access (OFDMA), peak-to-average 
T1 power ratio (PAPR). 



> 
■ ^ 

X 

H I. Introduction 

Orthogonal frequency division multiplexing (OFDM) is known as one of the best modulation 
schemes for high rate data transmission in wireless communication due to its robustness against 
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multipath fading, bandwidth efficiency, and simple implementation. However, due to high peak- 
to-average power ratio (PAPR) of the OFDM signals, the OFDM system requires expensive high 
power amplifiers having a large dynamic range. 

In the literatures, many PAPR reduction schemes for the OFDM systems have been proposed 
El-llEll- The simplest method for PAPR reduction of the OFDM signals is clipping [TJ. Clipping 
at the Nyquist sampling rate has been used for low complexity applications but suffers from peak 
regrowth after digital-to-analog (D/A) conversion. It is known that clipping of the oversampled 
OFDM signals reduces the peak regrowth after D/A conversion. But, it causes out-of-band 
radiation which has to be filtered. Furthermore, clipping of OFDM signals causes clipping noise 
which has sparsity in time domain. There are several schemes to mitigate this clipping noise in 
each of the two sampling rate cases ifTSl - lfTTll . The scheme in [15] shows good clipping noise 
cancellation performance, but it requires iterative maximum likelihood (ML) estimation for all 
tones with clipping and filtering at the receiver, which causes lots of computation. 

Compressed sensing (CS) is a sampling method that converts input signal in high dimension 
into the signal lying in the smaller dimension [|9l|-[[T2|. In general, it is not enough to recover 
an unknown signal using compressed observations in the reduced dimension. Nevertheless, 
if the input signal has sparsity, its reconstruction can be achieved at the receiver using CS 
reconstruction algorithm. In this context, the clipping noise at the receiver can be reconstructed 
by CS reconstruction algorithm. As the first work for this, tone reservation to reconstruct the 
clipping noise by CS is proposed in lfT3ll . where several tones are reserved at the transmitter 
before clipping and the receiver can reconstruct the clipping noise using CS reconstruction 
algorithm. However in this scheme, the reserved tones induce data rate loss and they should be 
demodulated at the receiver unlike the conventional tone reservation scheme. Additionally, this 
scheme shows poor bit error rate (BER) performance due to the small number of compressed 
observations (i.e., the small number of reserved tones) and vulnerability of CS reconstruction 
algorithm against noise. Motivated by lfT3l . [fT4l is proposed, where there is no data rate loss 
because pilot tones are exploited as compressed observations to reconstruct the clipping noise 
without reserved tones. But, it still shows poor BER performance for the same reasons as fTSl. 

Orthogonal frequency division multiple access (OFDMA) has severer PAPR problem due to the 
increased number of subcarriers. There are several reasons that many clipping noise cancellation 
schemes for OFDM systems such as the ones in [il31l and [|T6l cannot be directly applied to 
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the OFDMA systems. For example, the scheme in [[T5l requires iterative ML estimations for all 
tones together with clipping and filtering at the receiver, which is impractical for the OFDMA 
systems. For OFDMA systems, the clipping noise cancellation scheme is proposed in [fTTl . but it 
is not easy to be employed in the practical systems because the pilot tones should be estimated 
at the receiver but they are distorted by clipping. 

In this paper, we propose a new clipping noise cancellation method without data rate loss 
by using CS. The proposed scheme exploits compressed observations of the clipping noise 
underlying in the partial data tones to reconstruct the clipping noise without the reserved tones. 
In the proposed scheme, the number of compressed observations is adjustable in the sense that 
the optimal number of compressed observations can be selected according to the noise level. By 
doing so, we successfully overcomes the weakness of CS reconstruction against noise. Numerical 
analysis shows that the proposed method can mitigate the clipping noise well in both cases of 
clipping at Nyquist sampling rate and clipping and filtering at oversampling rate. Furthermore, 
the proposed scheme can also be applied easily to OFDMA systems with interleaved partitioning. 

This paper is organized as follows. The structure of OFDM system and clipping scheme 
are reviewed in Section II. In Section III, the proposed clipping noise cancellation scheme 
is presented. In Section IV, the application of the proposed algorithm to OFDMA systems is 
presented. Section V presents numerical analysis and the conclusion is given in Section VI. 

II. System Model 

Let X = (X(0),X(1), ...,X(A^ — 1))^ be an input symbol sequence. The continuous time 
baseband OFDM signal can be represented as 



where N is the number of subcarriers and T is the symbol duration. Let At^, = T/LN be a 
sampling interval, where L is oversampling factor. Then the discrete time OFDM signal sampled 
at time nAti can be expressed as 



An OFDM signal sequence with oversampling factor L can also be obtained by padding X with 
(L — 1)N zeros and processing the inverse discrete Fourier transform (IDFT). 




< t < T 



(1) 



n = 0, 1,...,LX- 1. 



(2) 
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The PAPR of the OFDM signal sequence XL{n) with oversampling factor L is defined as the 
ratio of the peak-to-average power of the signal as 

PAPR = ^^^o<n,.^ i |x.(n)P 
E[\xL[n)\^] 

where E\-] denotes the expected value. 

To reduce peak regrowth at the D/A conversion, clipping and filtering are performed on the 
OFDM signal sequence oversampled by a factor L at the transmitter. The clipping operation is 
given as 

xiin), \xL{n)\ < A 

XL{n) = I (4) 

\Aexp{jaig{xLin))}, |xL(n)| > A 
where A is the clipping threshold. Then the clipping ratio (CR) is defined as 

CR = 201og- [dB] (5) 
a 

where a = E[\xL{n)\'^]. 

It is shown in [|8l that the clipped signal xl^u) can be modeled in two different ways. The first 
one is an additive model where the clipped signal XL{n) is considered as the sum of the sampled 
signal XL{n) and the clipping noise CL{n), and the second one is an attenuated model where 
XL{n) is viewed as the sum of an attenuated component axL{n) and clipping noise component 
diin), that is, 

{Xhin) + cAn), for additive model 
(6) 
axi{n) + di{n), for attenuated model 

where n = 0,1, LN — 1 and the attenuation a is given as (HI 



a = 1 - e-^^ + ^^erfc(7) (7) 
and 7 = VlQCR/io. 

In order to remove the out-of-band radiation due to clipping operation, the clipped signal 
xi,{n) in time domain is converted back to the frequency domain by taking LA^-point discrete 
Fourier transform (DFT) and being filtered. Then we can have X{k) as 

A'(fc) = (8) 

aX{k) + -D(/c), for attenuated model 
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where A; = 0, 1, iV — 1. In [Sj, D{k) is assumed as a complex Gaussian random variable with 
zero mean. 

Let x{n) + c(n), n = 0,...,A^ — 1, denote the IDFT of ([S]) for the additive model, where 
x{n) is Nyquist sampled OFDM signal and c{n) is clipping noise. For the case of L = 1, we 
do not use filtering and c{n) is equal to ci{n). Then the clipping noise c{n) can be considered 
as K-sparse signal having K nonzero elements. We define the sparsity ratio as K/N . Now we 
can exploit this sparsity by CS method. If L is larger than 1, c{n) is not equal to CL{n) but c(n) 
can be considered as nearly /^-sparse signal, implying that the amplitudes of remaining N — K 
elements are close to zero, but not exactly zero. Even in this case, CS method to cancel the 
clipping noise can also be applied effectively. 

III. Clipping Noise Cancellation for OFDM Signals Using CS 
A. Formulation to CS Problem 

The received symbol using the additive model in ([8]) can be expressed in frequency domain 

as 

Y{k) = H{k){X{k) + C{k)) + Z{k), Q<k<N-l (9) 

where X{k) and Y{k) denote the input and received symbols, H{k) denotes the frequency 
domain channel response, and Z{k) denotes the additive white Gaussian noise (AWGN) with 
variance A^o in frequency domain. In matrix form, it can be rewritten as 

Y = H(X + C) + Z (10) 

where H = diag{H) and Y, X, C, Z are x 1 column vectors. 

We assume perfectly known channel response H and perfect synchronization. After channel 
equalization, its output is given as 

Y^""!^ = U-'^Y = X + C + U-^Z. (11) 

Note that in the attenuated model, Y^^'^^ can be similarly expressed as 

y(^'?) =aX + D + U^Z (12) 

where the clipping noise D is x 1 column vector. 
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Fig. 1. Block diagram of the proposed sciieme. 



For applying CS method, we need a compressed observation vector in the reduced dimension. 
Our suggestion here is to select a subset of components in F^"^"?), namely M out of N components 
of This can be done by multiplying M x N selection matrix Srr consisting of some M 



rows of the identity matrix I at to Y^^'^^ in ( 1 1 1. The discussion on S/jr will be made in the next 



subsection. Let C = Fc, where F is x DFT matrix. Then, we have 

Sr«F(^'') = SurFc + SurX + Srj,H-'Z. (13) 



If we subtract the estimation SrrX from ( |13| ), we have 

Y = SrrY^'''^ - SrrX = SrrFc + Srr{X - X) + SrrH-'Z 

= $c + Srr{X -X) + SrrH-'Z (14) 

V ' 

noise vector 

where the matrix $ = SrrF can be considered as M x measurement matrix in CS. As one 
can see in [fTOl . the measurement matrix for CS can be constructed by using the subset of rows 
in DFT matrix. Then the resulting equation ( [14] ) can be considered as CS problem, where the 
vector Y can be considered as M x 1 compressed observation vector, the clipping noise c as 
X 1 sparse signal vector, and the remaining vector Srr{X — X) + SrrYI^'^Z as M x 1 noise 
vector. By using CS reconstruction algorithm, we can reconstruct c as c from the compressed 
observation vector Y . Then fast Fourier transform (FFT) of c is subtracted from the equalized 
received symbol and then the final decision is made. Fig. [I] pictorially summarizes the 

proposed scheme. 
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B. Selection Matrix Srh 

As mentioned in the previous subsection, we suggest that the M x 1 compressed observation 
vector Y = SkrY^^'^^ - SrrX. And, SkrY^^'^^ in Y is constructed from Y^'"'^ by selecting those 
components of F^'^'J) ^j^^t are thought to be more reliable than others. We designate a specific 
area in the received signal space as reliable region (RR) and define the index set ICrr as the 
set of component indices of the received signals that lie in the reliable region. For quaternary 
phase shift keying (QPSK) and 16-quadrature amplitude modulation (QAM), RR is designated 
as shaded area in Fig. |2] For other modulations, RR can be set in analogous way. In Fig. [2[ the 
distance between adjacent signal points is assumed to be 2. 




Fig. 2. Reliable regions for QPSK and 16-QAM modulations. 



Prior to explain the index set ICrr in detail, we compare the additive model and the attenuated 
model in ^ from the viewpoint of signal-to-clipping-noise ratio (SCNR). In (HI, D{k) in 
attenuated model is uncorrelated clipping noise from X{k). And, C{k) in additive model is 
easily derived from ([8]) as 

C{k) = {a-l)X{k) + D{k). (15) 



Additionally to D{k), C{k) in ( 15 1 also has a term {a — l)X{k) correlated from X{k). And thus 
it implies that the SCNR of the attenuated model is larger than SCNR of the additive model. 
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that is, 

E{\aX{m E[\X{m 

E[\D{m - E{\c{m- ^ ^ 

For the reason as above, we describe the index set ICrr in detail using the attenuated model. 
Assume that M components of {l/a)Y^^'^^ fall into RR. Then the index set ICrr can be expressed 
as 

ICrr = {k^ : -Y^''^\km) G RR, < m < M - 1}. (17) 

a 

The M X N selection matrix Srr can be obtained from the identity matrix of order by 
selecting M rows corresponding to the index set ICrr, where M is cardinality of the set ICrr. 

In our proposed scheme, for the determined Srr, we estimate SrrX. Clearly, we can also 
conclude that SrrX can be more efficiently estimated by using attenuated model. That is, in 
order to estimate SrrX efficiently, we decide (l/a)Sij/jF*^'^'^^ as SrrX by ML estimation. 
Intuitionally, at the receiver, when we decide {l/a)Y^'^'^\km) in RR as X{km), its decision error 
probability is lower than that when we decide (l/a)Y^^'^\k ^ ICrr) in the outside region of RR 
as X(A; ^ ICrr). Therefore, if we consider the signals only in the RR to reconstruct the clipping 
noise c, decision error of SrrX can be reduced as 5 increases in Fig. |2| 

C. Decision Error Probability of SrrX Over an AWGN Channel 



For an AWGN channel, after scaling by 1/a, (12) can be rewritten as 



-Y{k) = -Y'^^'^\k)=X{k) + -W{k), fc = 0,l,...,A^-l (18) 

a a a 

where W{k) = D{k) + Z(k). Since D(k) is assumed to be complex Gaussian random variable 



with zero mean, {l/a)W{k) in (18) can be characterized as a zero mean complex Gaussian 
random variable with variance E[\{l/a)W{k)\'^] = {l/a'^){E[\D{k)\'^] + Nq), statistically inde- 
pendent of the input symbol [8|. We only consider the case of clipping at Nyquist sampling rate. 
For clipping at Nyquist sampling rate, all the data tones experience the clipping noise D{k) with 
the same variance as 

E[\D{k)\'] = (1 - e-^' - a')E[\X{k)\']. (19) 
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And thus the variance of {l/a)W{k), denoted by Nq = E[\{1 / a)W {k)\'^] is the same for each 
data tone. The decision error probability of SrrX for V^-QAM of even V is expressed as 



a 



Pv2 (X{k) ^ X{k) I -Y{k) E RR) = I - Py2 [X{k) = X{k) \ -Y{k) G RR 



a 



1 



Py2(X(A;) = X{k) n lY{k) e RR) 
Pv^&{k) ERR) ' 



(20) 



In (|20|), Pv2{{l/a)Y{k) G RR), the probability that {l/a)Y{k) belongs to RR for V^.qam 
is rewritten as 



a 



Pv2[-Y{k) E RR) = Pv[-Y{k) E RR 



1 



a 



(21) 



where PY({l/a)Y(k) E RR) is the probability that inphase (or quadrature) component of 
{l/a)Y{k) belongs to RR for V-pulse amplitude modulation (R^M). 
Then, PY((l/a)Y(k) E RR) is expressed as 



Pv{ -Y{k) E RR 



(y-2)/2 

v=0 
A-V-5 



V 

+ 

+ ... + 

o (V'-2)/2 



2~V-S 



-{x-{2v+i))yN;, 



1 



2-V+5 a/ttA^ 
2+V-5 



-{x-{2v+l))yN[, _^ 



y-V-5 



4-V+5 a/ VtA^o 



-{x-i2v+l)y/% 



v=0 
V-3 I 



2t; - 1 + 1/ + 5 



+ g 



2t;-3 + + 



-{x-{2v+\)f 1% 



■2v + 2v' + 1-V + 5 



Q\ 



-2t; + 2t;' + 3 - r - 5 



(22) 



where (5(2) = /. 



The probability Py2(X(/c) = X(A;) n {l/a)Y{k) E RR) in (|20|) is equivalent to the correct 
decision probability for V^-QAM when distance between adjacent signal points is 2 — 25 and it 
is given in [|20ll as 



Pv2(x{k) = X{k) n -Y{k) E RR 



a 



1-5 



(23) 



Consequently, the decision error probability of SrrX is obtained by substituting (22) and (23) 
into (I20I). 
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D. Determination of 6 

Accuracy of reconstruction via CS in ([14]) is affected by not only the noise but also the 
selection of compressed observations. Certainly, both the noise and the compressed observations 
vary according to 5 in our scheme. The influence of 5 on the performance of CS reconstruction 
can be explained in two different views. 

Firstly, the amount of estimation errors is dependent on the value of 5. If 5 is set to be 
too small, then estimation error will occur frequently, which results in increasing the energy of 



|S/jij(X — X)| in the noise vector in (14). Accordingly, the performance of CS reconstruction 
will be degraded and thus we can say that large value of 5 is preferable in this point of view. 

Secondly, the number of compressed observations which is the row size M of $ = SrrF is 
dependent on 5. As 5 decreases, the RR becomes wider, which implies that M becomes larger 



as 



E[M] = N ■Pv2 {^y{k) G RF^j . (24) 

In our proposed scheme, for clipping at Nyquist sampling rate over an AWGN channel, the 
probability that the component k is in ICrr is identical for all k. Accordingly, $ can be considered 
as random selection of M rows from A^-point DFT matrix. In [11 J, for randomly selected 
compressed observations in frequency domain, the quality of the reconstruction is upper bounded 
as 

\\c-ch<0{l)-R-{M/{\ogNf)-\ (25) 

where r and R are related to the sparsity of the clipping noise c. And thus, for fixed N and 
c, large M improves the reconstruction performance. Therefore, in our proposed scheme, small 
value of 5 is preferred in this point of view. 

Consequently, we have to consider above two aspects for determining 5. We perform the 
massive simulations to find an optimal value of 5 which maximizes the performance of CS 
reconstruction. 



E. Reconstruction Algorithms in CS 

A number of reconstruction algorithms have been introduced in many literature of the CS 
ll9]l- [[TT]| . We adopt a sparse approximation algorithm called orthogonal matching pursuit (OMP) 
to handle the signal recovery problem [TS]. The main reason of selecting OMP algorithm is 
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not its performance but its ease of implementation and speed. Other reconstruction algorithms 
such as basis pursuit (BP) can show better reconstruction performance, while they require high 
computational complexity. 

Since there is no information about sparsity K of clipping noise c at the receiver, we need to set 
the iteration number of OMP sufficiently large. As OMP algorithm runs iteratively, it estimates 
the nonzero clipping noise in descending order of amplitude. The clipping noise having small 
amplitude below a certain threshold can be negligible. Therefore, at each iteration of OMP, when 
the recent reconstructed clipping noise component is below the properly predetermined threshold 
or the maximum number of iterations is reached, we can stop the OMP algorithm to reduce the 
unnecessary computational complexity. 

In this paper, we investigate into the case of clipping and filtering (L > 1) and the case of 
clipping at Nyquist sampling rate (L = 1). Additionally, the K largest peaks reduction scheme 
similar to the peak annihilation in [13] is also simulated. That is, we scale down exactly K 
components having large amplitudes in the OFDM signal sequence x at Nyquist sampling rate. 
These components are scaled down to have 80% of (K + l)-th largest OFDM signal's amplitude 
to mitigate peak regrowth after D/A conversion. In this case, the OMP needs to be done K 
times. 

IV. Clipping Noise Cancellation for Interleaved OFDMA Downlink Systems 

Different from OFDM systems, in the downlink of OFDMA systems, each user only has 
knowledge of its own data tones. Moreover, users in OFDMA system can use different modulation 
schemes in general. Since it is not always possible for a particular user to demodulate the data 
tones of the other users, the estimation SjirX of the data tones in the RR could be impossible at 
the receiver. Hence the proposed scheme needs to be modified in order to be applied to OFDMA 
systems. 

In a 2^-user OFDMA downlink system with N = 2^ subcarriers, all the available subcarriers 
can be partitioned into 2^ interleaved subcarrier groups. For each subcarrier group assigned to 
each user, small FFT block can be extracted from the point FFT block. The main idea of this 
section is exploiting this structure. 

Fig. |3] shows the structure of decimation-in-frequency FFT of the clipping noise c when = 8 
and U = 1. As one can see in Fig. [3| subcarriers are assigned to user 1 and user 2 by interleaved 
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pattern. In general, A^-point FFT consists of b stages and then this can also be represented by 
the matrix form as 

F = OT,T;,_i...T2Ti (26) 

where Tj, I < i < b, denotes the i-th stage matrix of FFT process and O is an index ordering 
matrix. For instance, 8-point FFT block can be represented by the matrix form OT3T2T1 as 
shown in Fig. [s] Similarly, for A^/2^-point FFT, we also define the matrices O' and T'-,, 1 < 
i' < b — U, as the index ordering matrix and i'-th stage matrix of FFT process, respectively. 
Let Iat be the identity matrix of order and ej be the j-th row of Iat. Let and be 
X interleaved and adjacent partitioning matrices for the u-th user defined as P^ = 
[e^e^[/^^...e^_2(7^„]^ and P;^ = {D'y^P^O. For instance, in Fig. [3| for user 1, P{ selects the 
subset of components in the clipping noise C (after O) denoted as "user 1" (i.e., C(0), C(2), C(4), C(6)). 
And, P^ selects the subset of components between Ti and T2 denoted as "user 1" (i.e., c'(0), 
c'(l), c'(2), c'(3)). 

In order to extract the A^/2^-point FFT block for the u-th user from the -point FFT block 



for OFDMA signals, we multiply (11 1 by the u-th user interleaved partitioning matrix P^. And 



December 9, 2011 



DRAFT 



IEEE TRANSACTIONS ON SIGNAL PROCESSING 



13 



using ( [26] ), (fTTI) can be modified as 

p/y(eg) ^ pi^x + P^(OTfe...Ti)c + PiH^Z 

= PiX + Pi{OTk...Tu+i)iTu...T,c) + PiH'Z. 



(27) 



Clearly, P^O can be modified as O'P;^ and it is easy to check that P^{Tb...T 



u+i) 



{T[_jj...T[)P^. Then (27) can be rewritten as 



(28) 



p/y{e,) ^ p/^ ^ 0'P^{T,...Tu+i)iTu...T^c) + PiWZ 

= PiX + o'(tU-t;)p;!(Tc/...Tic) + Piu-'z 

= PiX + F'c' + PiH^Z 

where F = 0'{T'f,_u...T[) is x DFT matrix and c' = P^Tu-.-Tic. 

Since x matrix P^T[/...Ti has only one nonzero elements per a column, if c is 

K sparse, then c' = P^Tu---Tic is also sparse signal having the maximum K sparsity but the 
reduced dimension from N to Therefore, the sparsity ratio of clipping noise c' in the 



OFDMA system effectively increases compared to c in the OFDM system in (14). And, it may 
somewhat degrade reconstruction performance. 



Now, as we have previously done in (13) and (14) for OFDM systems, we can formulate (28) 
to the CS problem as 



Y' = S'^^P^F(^«) - S'^^P^X = S'^^F'c' + S'^nKiX -X) + S'^^P^R-^Z 

= $V + S'^^P^X -X) + S'^^P^H-^Z 



(29) 



noise vector 



where S'j^j^ and 



S'^^F' are selection matrix and a measurement matrix for OFDMA 



systems, respectively. In (29), it is possible to estimate S'j^j^P^X because it requires only the 
data tones of ?i-th user. In the proposed noise cancellation scheme for OFDMA systems, from 
the compressed observation vector Y', we reconstruct the clipping noise c', which is associated 
with only the data tones of u-th user, instead of reconstructing the whole clipping noise c. For 
instance, in Fig. [3| clipping noise c'(0), c'(l), c'(2), c'(3) is reconstructed for user 1 instead of 
clipping noise c(0), c(l), c(2), c(7). 
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V. Numerical Analysis 

In this section, we evaluate the BER performance of the proposed clipping noise cancellation 
scheme in the AWGN channel. We simulate both the case of clipping at Nyquist sampling 
rate with L = 1 and the case of clipping and filtering at the oversampling rate with L = 4. 
Additionally, the K largest peaks reduction scheme is also simulated. The distance between 
adjacent signal points is 2. 




SNR [dB] 



Fig. 4. BER performance of the proposed scheme and the scheme in 1131 with K largest peaks reduction when A'^ = 64 and 
QPSK is used. 

Fig. |4] shows the BER performance of the proposed scheme using K largest peaks reduction. 
There is a large benefit to use the proposed scheme compared to no clipping noise cancellation 
case. As K increases, the BER performance decreases and for > 8, the BER performance is 
not closed to that of the original OFDM. Thus, we use the maximum iteration number of OMP 
as 0.125A^ for unknown K cases. In Fig. |4} clipping noise cancellation scheme in [[T3l is also 
given when 29 tones are positioned by a (59, 29, 14) difference set [fT9l and reserved. It uses the 
sufficient number of reserved tones (i.e., the sufficient number of compressed observations) for 
reconstructing the clipping noise with K = A in the absence of noise. But, the scheme in llT3ll 
shows the poor BER performance due to weakness of CS reconstruction against noise. Although 
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the clipping noise cancellation scheme in [[T4]| is not shown in Fig. |4} it may show worse BER 
performance than that in [[T3l. The reason is that, considering the practical system, the number 
of pilot tones (i.e., the number of compressed observations) is smaller than 29 when = 64. 
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Fig. 5. BER performance of the proposed scheme with chpping at Nyquist samphng rate when CR = OdB, — 64, and 
QPSK is used. 



Fig. |5] shows the BER performance of the proposed scheme using clipping at Nyquist sampling 
rate for various 6 when = 64 and QPSK is used. We set the maximum iteration number for 
OMP as 8. It shows that 5 = 0.4 ~ 0.5 gives the best BER performance. As we mentioned, 
smaller or larger value than S = 0.4 ~ 0.5 causes frequent decision error or lack of compressed 
observations. Thus it induces the degradation of BER performance. In Fig. [5| since CR= OdB 
is small value, K nonzero components of the clipping noise c is larger than 8. Therefore, 
OMP cannot reconstruct the remaining nonzero components in the clipping noise c and BER 
performance is not closed to that of the original OFDM. Additionally, for S = 0.0, we do not 
need to use CS reconstruction algorithm because M = N. Using largest 8 components in IDFT 
of the compressed observation vector for reconstructed clipping noise c gives us almost the same 
BER performance as using CS reconstruction. 
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Fig. |6] shows the BER performance when clipping and filtering are performed (L = 4) and we 
reconstruct the clipping noise c by using OMP for several S, N, and modulation schemes. We 
set the maximum iteration number for OMP as 0.125A^. Similar to Fig. [5| there is a 5 which can 
achieve the best BER performance. Since clipping noise c for L > 1 has nearly sparsity, it shows 
the limited performance compared to clipping at Nyquist sampling rate. Fig. |6] also shows the 
BER performance of the proposed scheme for large N. For large N, the BER performance is 
almost the same as that of small A^. And, 5 for the best BER performance is the same. Therefore, 
we can conclude that 5 for the best BER performance is mainly affected by CR and modulation 
order. 
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Fig. 7. BER performance of the proposed scheme applied to OFDM A system for 4 users with clipping at Nyquist sampling 
rate when CR = 2dB, = 512, and QPSK is used. 



Fig. |7] shows the BER performance when the proposed scheme is applied to interleaved 
OFDMA systems. Clipping at Nyquist sampling rate is used for = 512 and 2^ = 4 users. And 
we set the maximum iteration number for OMP as 16 and 6 = 0.2 is used. We consider only 
the 128 data tones corresponding to user 1 and also construct only the clipping noise for user 
1. Therefore, its computational complexity is reduced compared to OFDM system for A^ = 512 
because the size of measurement matrix decreased. However, sparsity ratio increases by the 
maximum 4 times and it causes poor reconstruction performance. 

VI. Conclusion 

In this paper, we proposed the new clipping noise cancellation scheme in OFDM using CS. 
In the proposed scheme, the data tones are partially exploited to reconstruct the clipping noise 
instead of the whole data tones. By introducing RR, we can select these partial data tones 
and numerical analysis shows that the proper value of 6 needs to be chosen. Additionally, the 
proposed scheme can also be applied easily to OFDMA systems because of the FFT structure. 
Using the proposed clipping noise cancellation scheme for the OFDM and OFDMA systems, 
the BER performance can be improved compared to the conventional schemes. 
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